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Analysis  of  Free  Electron  Laser  Performance  Utilizing 
the  National  Bureau  of  Standards'  CW  Micro tron 

Introduction 


An  FEL  facility  for  applications,  primarily  in  biomedical  and  material  science  research 
as  well  as  for  basic  physics  and  chemistry,  is  to  be  situated  at  the  National  Bureau  of 
Standards1,2.  A  CW  185  MeV  racetrack  microtron  (RTM)3  is  under  construction.  The 
NBS  Accelerator  Laboratory  consists  of  a  series  of  interconnected,  individually  shielded, 
underground  halls.  The  updated  layout  is  indicated  in  Fig.  1.  The  FEL  is  expected  to  be 
operational  by  1990. 

The  major  limitation  of  an  RTM  as  an  FEL  driver  is  that  its  peak  current  capability 
is  lower  than  electron  linacs  which  operate  in  the  same  energy  range.  However,  the  RTM 
is  superior  to  pulsed  linacs  in  energy  spread  and  emittance.  The  RTM  is  comparable  to  a 
storage  ring  in  terms  of  beam  emittance  and  energy  spread,  but  there  is  no  restriction  on 
insertion  length  or  “stay  clear”  aperture.  The  beam  energy  can  be  varied  continuously  over 
a  wide  range  without  significant  loss  of  performance.  In  addition,  microtrons  are  compact 
and  energy  efficient.  Because  of  the  CW  nature  of  the  RTM,  the  generation  of  coherent 
photons  is  not  hindered  by  a  finite  macropulse  length. 

The  original  design  parameters  of  the  NBS  RTM  are  given  in  Refs.  1-3.  The  design 
calculations  indicate  a  longitudinal  emittance  <  30  keV-degrees  and  a  normalized  trans¬ 
verse  emittance1,2  e,v  <  10  mm-mrad.  Based  on  recent  measurements  of  the  performance 
of  the  5  MeV  injector  linac,  the  actual  values  of  both  the  longitudinal  and  transverse 
emittance  are  expected  to  be  smaller  than  the  design  values.  The  injector  system  must 
be  upgraded  to  provide  a  peak  current  of  >  2  A  in  3.5  psec  micropulses,  giving  electron 
pulse  length  le\>  ~  0.1  cm.  In  order  to  keep  the  average  electron  beam  power  within  the 
capability  of  the  existing  RF  power  system,  the  new  injector  will  fill  only  a  small  fraction  of 
the  RF  buckets  (e.g.,  1/24,  1/120  depending  on  electron  beam  energy).  We  are  proceeding 
with  a  design  of  a  photocathode  injector  system  lor  this  upgrade. 

1-P  Fret1  electron  Laser  Analysis 

A  first  order  evaluation  of  the  I'LL  p*-i  fin  im.um  •  •  "I  He  NM'-i  HIM  •  an  be  obtained 
from  the  1-D  small  signal  low  gain  formula4.  The  results  indicate  that  sufficient  gain  can 
Manuscript  Approved  October  30,  1987. 
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be  obtained  at  fundamental  wavelengths  in  the  range  from  10  /im  >  A  >  0.25  /xm.  The 
formula  for  the  electric  field  amplitude  gain  G  in  the  small  signal,  low  gain  regime,  can  be 
written  as 

°  =  T  tfK'N'r  (  — y. 

(TR  /  4  j"  du  \  v  ) 

where  N  is  the  number  of  wiggler  periods,  -ya  is  the  initial  relativistic  gamma  factor, 
=  7r r2  is  the  cross-sectional  area  of  the  radiation,  r0  is  the  minimum  1/e  radius  of 
the  Gaussian  radiation  field  amplitude,  IA  =  17  x  103  A,  I  is  the  current  in  Amperes, 
K  =  (\e\Bw\w/2wm0c2)R1\fs  is  the  wiggler  parameter,  Bw  is  the  magnetic  field  in  the 
wiggler,  Au,  is  the  wavelength  of  the  wiggler,  F\  =  J0(b )  —  J i(b)  for  a  linearly  polarized 
wiggler,  b  =  A'2 /2(  1 +  K 2 ),  v  —  —  N\{uj  —  u)0)j2c  is  the  normalized  frequency  mismatch,  and 
u>0  ~  2y2c(  27r/Au, )/( 1 -I- A'2 )  is  the  resonant  angular  frequency.  The  function  d/ di/{s\n  u /  u)2 
has  a  maximum  value  of  0.54  when  v  =  —1.3. 

The  power  gain  can  be  obtained  by 

Gp  =  (1  +  G)2  -  1. 

In  the  low  gain  regime,  Gp  ~  26’.  The  FEL  will  oscillate  when  the  power  gain  is  greater 
than  the  losses  per  pass  in  the  reson-ator.  The  1-D  gain  formula  is  only  a  rough  estimate. 
It  is  sensitive  to  the  choice  of  filling  factor. 

The  conceptual  design  consists  of  a  linearly  polarized  wiggler  with  a  period  of  Au,  =  2.8 
cm,  and  a  nominal  magnetic  field  amplitude  of  Bwo  —  5400  G.  This  can  be  constructed 
with  a  hybrid  wiggler  design  with  the  gap  separating  the  wiggler  poles  of  g  =  1.0  cm. 
A  wiggler  can  be  constructed  conceptually  in  more  than  one  section,  such  that  a  wiggler 
of  shorter  length  can  also  be  available.  A  shorter  wiggler  and  a  corresponding  vacuum 
chamber  may  be  necessary  for  long  wavelength  operation. 

Figure  2  is  a  plot  of  the  1-D  maximum  small  signal  power  gain  versus  wavelength, 
assuming  a  conservative  peak  current  of  2  A.  The  open  circles  !>>)  are  obtained  with 
electron  beam  energies  of  25.  50.  i  5.  [25.  and  I  <  5  Me\  .  Tin’  solid  curves  are  obtained  tor 
the  same  electron  beam  energies,  but  varying  He  \  iggbu  ampin  mb-  Icmi  tbb/d,,.,,  to 
the  magnetic  field  in  the  wiggler  is  to  be  changed  by  varying  the  gap  between  the  poles 


from  1.4  cm  to  1.0  cm.  As  the  magnetic  field  decreases,  the  wavelength  of  the  radiation 
decreAses,  and  the  gain  is  reduced. 

3-D  Effects  on  the  Gain 

Since  FELs  are  not  actually  1-D,  3-D  effects  will  change  the  gain.  Some  of  the  3- 
D  effects  that  we  will  examine  in  this  paper  are  finite  transverse  emittance,  radiation 
diffraction  and  refraction,  and  some  effects  associated  with  finite-length  electron  pulses. 

We  will  assess  these  three-dimensional  effects  using  a  fully  3-D  self-consistent  computer 
code,  SHERA,  developed  at  the  Naval  Research  Laboratory.  The  formulation  of  the  wave 
equation  is  based  on  the  source  dependent  expansion5  of  the  radiation  field,  and  the  electron 
dynamics6  are  evaluated  self-consistently.  We  assume  a  waterbag  distribution  in  the  4-D 
transverse  emittance  space,  which  leads  to  a  parabolic  profile  for  the  electron  beam  density. 
Since  the  energy  spread  of  the  NBS  RTM  is  very  small,  it  will  not  be  considered;  and  we 
will  also  not  treat  the  effects  of  pulse  slippage  on  the  gain.  The  radiations  are  taken  to 
have  a  Rayleigh  length  of  175  cm  with  the  minimum  radiation  waist  located  at  the  center 
of  the  wiggler.  Results  for  two  different  operating  regimes  will  be  presented. 

The  effect  of  the  emittance  on  the  performance  of  the  FEL  will  be  more  important 
for  short  wavelength  operations.  Thus,  our  first  example  will  be  for  A  =  0.23  /am  with 
70  =  350.  The  pulse  slippage  distance,  N A  =  0.003  cm,  is  much  shorter  than  the  electron 
pulse  length,  .  The  minimum  1/e  radiation  field  amplitude  waist  is  ra  =  3.57  x  10-2 
cm.  Plots  of  the  power  gain,  Gp,  versus  the  normalized  frequency  mismatch,  u,  are  shown 
in  Figure  3.  Curve  (a)  gives  the  1-D  estimate  of  the  gain.  Curves  (b),  (c)  and  (d)  are  the 
gains  calculated  from  the  computer  code  for  normalized  transverse  emittance  of  e^r  =5, 
10  and  20  mm  inrad,  respectively.  The  radii  of  the  electron  beams  were  determined  by 
properly  matching  the  beam  into  the  wiggler,  i.e.,  the  radii  of  the  beams  inside  the  wiggler 
is  uniform.  The  matched  beam  radius  condition  is 


where  I\j  =  \J2tc  K  /  \  is  the  betatron  vvav<nuinl»r  (or  with  parabolic  pole  laces', 

where  the  focusing  in  both  transverse  directions  is  equal.  The  matched  edge  radius  of  the 
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independent  of  the  beam  energy.  The  matched  edge  radii  of  the  electron  beam  are  reb  = 
1.77  x  1 0 —  2 ,  2.50  x  10-2  and  3.54  x  10-2  cm  for  normalized  edge  transverse  emittances  of 
ey  —  5,  10  and  20  mm  mrad,  respectively.  If  the  emittance  becomes  larger  than  20  mm 
mrad,  the  radius  of  the  electron  beam  will  become  larger  than  the  radiation  spot  size,  and 
the  gain  will  be  substantially  reduced. 

The  effect  of  finite  emittance  on  the  gain  is  negligible  for  A  =  1.25  pr n  with  =  150. 
The  pulse  slippage  distance  in  this  case  is  0.016  cm,  and  it  is  still  unimportant.  Figure  4 
shows  plots  of  the  power  gain,  Gp,  versus  the  normalized  frequency  mismatch,  i/,  similar 
to  Fig.  3.  Again,  the  curve  (a)  gives  the  1-D  estimate  of  gain.  Curves  (b),  (c)  and  (d) 
are  the  gains  calculated  from  the  computer  code  for  normalized  emittance  of  epj  =5,  10 
and  20  mm  mrad,  respectively.  Since  the  wavelength  is  longer,  the  minimum  l/e  radiation 
field  amplitude  waist  becomes  ra  =  8.3  x  10~2  cm,  and  the  electron  beam  radii  are  much 
smaller  than  the  radiation  waist.  The  gain  at  A  =  1.25  pm  is  insensitive  to  the  design 
value  of  the  finite  transverse  emittance. 

Figures  3  and  4  also  show  a  shift  of  the  zero  crossing  of  the  gain  curves  obtained  from 
3-D  simulation.  This  shift  comes  from  the  change  in  the  phase  of  the  diffracting  radiation 
field.  It  has  no  real  important  effect  on  the  oscillation  criteria  for  the  examples  under 
consideration. 

Figure  5  plots  the  maximum  3-D  power  gain  versus  wavelength  with  a  peak  current  of 
2  A,  for  normalized  emittances  of  5,  10,  and  20  mm  mrad.  Each  curve  is  obtained  for  the 
identified  electron  beam  energy,  but  varying  the  magnetic  wiggler  amplitude  from  Bwo  to 
0.6 Bwo,  where  the  longer  wavelength  correspond  to  the  larger  magnetic  field.  Normalized 
emittance  is  very  good  in  the  long  wavelength  operating  regime.  In  the  shortest  wavelength 
operating  regime,  the  normalized  emittance  larger  than  5  mm  mrad  should  be  avoided. 


Since  the  current  is  a  function  ol  axial  p<>  Hem  m  a  lim'<-  fi  ngtli  electron  pulse,  and 
pulse  slippage  is  unimportant  ,  the  local  gain  is  a  function  of  the  local  current  in  the  electron 


pulse.  For  the  first  example  at  A  =  0.23  /xm  with  normalized  transverse  edge  emittance 
=  10  mm  mrad,  the  simulations  indicate  that  the  gain  is  proportional  to  the  local 
current,  consistent  with  the  1-D  formula.  For  the  second  example  at  A  =  1.25  fim,  the 
gain  increases  faster  than  the  linear  power  of  the  current.  Figure  6  is  a  plot  of  normalized 
power  gain,  i.e.,  power  gain  from  simulation  divided  by  the  maximum  1-D  power  gain, 
versus  the  frequency  mismatch  at  A  =  1.25  with  normalized  transverse  edge  emittance 
tjV  —  10  mm  mrad  for  three  different  currents:  (a)  /  =  4.0  A,  (b)  I  =  2.0  A  and  (c)  /  =  0.5 
A.  We  find  that  the  normalized  gain  increases  as  current  increases.  This  can  be  explained 
by  the  self-focusing  or  guiding  phenomenon5,8-11  of  the  FEL.  This  is  most  easily  observed 
in  the  plots  of  the  normalized  1/e  Gaussian  radiation  field  amplitude  radius,  shown  in  Fig. 
7.  For  I  =  0.5  A,  the  radiation  radius  behaves  like  a  free  space  resonator  radiation  field, 
curve  (-  -  -).  For  I  =  2  A,  the  radiation  radius  is  less  than  the  free  space  radius  at  the  end 

of  the  wiggler  as  self-focusing  begins  to  show,  curve  ( - ).  If  the  current  can  be  increased 

to  4  A,  the  radiation  becomes  even  more  focused,  curve  ( . ).  The  reason  that  self-focusing 

is  evident  at  such  low  current  is  that  the  emittance  is  very  good  and  current  density  is  high 
throughout  the  interaction  region,  i.e.,  high  beam  brightness  Bn  =  2l/(-K2e2N)  >  4  x  109 
A/nr/rad2,  where  I  >  2  A  and  edge  emittance  c,\-  =  10  mm  mrad. 

Conclusions 

The  3-D  self-consistent  simulation  results  from  the  computer  code  SHERA  indicate 
that  the  design  value  of  the  transverse  emittance  is  very  good,  so  that  it  does  not  degrade 
the  FEL  performance  for  intermediate  to  long  wavelengths.  For  the  shortest  wavelength 
under  consideration,  emittance  larger  than  5  mm  mrad  should  be  avoided.  Due  to  the 
good  emittance,  the  current  density  is  high  enough  that  focusing,  or  guiding,  begins  to 
manifest  itself  for  wavelengths  2.0  (mi. 
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Fig.  3.  Power  gain,  Gp ,  versus  frequency  mismatch,  v,  at  A  =  0.23  fim  with  = 
350.  Curve  (a)  is  based  on  1-D  gain  formula.  Curves  (b),  (c)  and  (d)  are  obtained  from 
simulations  with  normalized  transverse  edge  emittances  of  e/v  =  5,  10  and  20  mm  mrad, 
respectively. 
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Fig.  4.  Power  gain,  G'p,  verus  frequency  mismatch,  t/,  at  A  =  1.25  ptm  with  7„  = 
150.  Curve  (a)  is  based  on  1-D  gain  formula.  Curves  (b),  (c)  and  (d)  are  obtained  from 
simulations  with  normalized  transverse  edge  emittances  of  tjv  =  5,  10  and  20  mm  mrad, 
respectively. 


Radiation  Wavelength  X  (^m) 

Fig.  5.  Power  gain  versus  wavelength  based  on  a  fully  3-D  self-consistent  simulation  by 
varying  energy  and  einittance  of  the  electron  beam,  and  the  magnetic  field  of  the  wiggler. 


Fig.  6.  Normalized  gain  versus  frequency  mismatch  for  A  =  1.26  n m  and  edge 
emittance  of  €/v  =  10  mm  mrad.  Curves  (a),  (b)  and  (c)  correspond  to  results  obtained 
with  currents  of  /  =  4.0,  2.0,  and  0.5  A,  respectively. 


Axial  Distance,  z 

Fig.  7.  Normalized  1/e  Gaussian  radiation  field  amplitude  radius  as  a  function  of 
distance  z  in  the  wiggler  with  A  =  1.25  /mi  and  edge  emittance  e^v  =  10  mm  mrad  for 
three  different  currents:  ( - )  /  =  0.5  A.  { - )  I  -  2  A  and  ( . )  I  —  1  A. 
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